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Abstract 

Formulas for matrix determinants, algebraic adjunctions, characteristic polynomial co- 
efficients, components of eigenvectors are obtained in the form of signless sums of matrix 
elements products taking by special graphs. Signless formulas are very important for singu- 
lar and stochastic problems. They are also useful for spectral analysis of large very sparse 
matrices. 

1 Introduction 

Graph theory is a natural instrument for matrix determinants and combined features 
calculating. It is clear that determinants can be expressed in form of alternating sum of 
matrix elements products taking by some graphs corresponding to index permutation, 
each term sign is determined by permutation even. However, the alternation itself is 
an essential obstruction both in numerical calculating and in analytic research of spec- 
trum properties. It concerns for example cases when matrix elements depend on small 
(large) parameter and possess on it different orders. Thus under singular perturbed 
equations research |^ (in particular such as Fokker-Plank equations) necessity appears 
for spectral analysis of big order matrix having exponentially small elements |]^, |^. The 
same situation appear under consideration of Markov's chains connected with diffusion 
process where the questions of stochastic continuity, possible subprocess kinds and 
their structure properties run into the necessity of spectral analysis just such kind of 
matrixes. 

Our aim is to obtain the formulas for matrix determinants, algebraic adjunctions, 
characteristic polynomial coefficients, components of eigen-vectors in the form of sign- 
less sums of matrix elements products taking by special graphs. It turns out that one 
can get such formulas in terms of so called "tree" structure of some graph corresponding 
to matrix. 

The first step in this direction was made by Kirchhoff who computed the num- 
ber of connected subgraphs containing all vertices and containing no circuits (spanning 
trees). This number turned out to be equal to cofactor of any element of so called 
conductivity or Kirchhoff matrix of non-directed graph. Later this theorem was gener- 
alized for all coefficients of characteristic polynomial of this matrix, and also for cases 
of directed graphs where every arc (ordered pairs of vertices) possesses some quantity 
called weight Here we get the same type formulas not for special matrices (as 
Kirchhoff type) but for arbitrary ones. 



2 Main definitions and designations 

Unification of designations and even terminology proper is not complete yet in graph 
theory. So firstly we adduce the necessary definitions and notations. 
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Let G be digraph (directed graph). We use VG and AG to denote the set 
of vertices and arcs of G . The subgraph if of G is called factor if VH = VG. The 
outdegree (indegree) of the vertex i (the number of arcs going out of (into) i ) we 
denote d'^{i) {d~{i)). 

If G is digraph in which every arc has its own weight Qij (weighted adjacencies 
digraph), corresponding matrix G = {gij}i^j=i is called generalized adjacency matrix 
(the element gij = if there is no arc in G). 

The sequence of following each other arcs along their orientations is called a way if 
all vertices besides possibly the uttermost ones are different. The way connecting the 
vertices m and n we denote m ■ n. The cyclic way is called dicircuit. Linear digraph is 
digraph every vertex of which has unit in- and out-degree. So it consists of dicircuits. 

We associate with any weighed adjacencies digraph G the quantity tcg by the rule 

ttg = n 9tj 

which is naturally to call by productivity. 

Later on forests are the main graph theory object we use. As known there are two 
forest kinds in digraph situation. Here we call by forest digraph without dicircuits in 
which every vertex outdegree is equal to or 1 {d~^{i) = 0, 1). The only vertex of tree 
(component of forest) having zero outdegree {d'^{i) = 0) we call root. 

Let Af be finite set and W is some subset of A/". By J^^, k < \Af \ - \ W\ we denote 
the set of forests F obeying the following conditions 

1) VF = Af; 

2) F consists of exactly k + |VV| trees; 

3) The set of roots of F contain W ClS db subset. 

Suppose also JF'^ = jFg , jFyy = JF-J^. Note that the set is empty set and JF'-'^I consists 
of the only empty forest having only roots and no arcs. 

If some additional condition z is put on forests from JF^ we denote such set of 
forests as Tif{z). Other necessary utilized notations we sign as necessary and also we 
use sometimes the term "graph" in wide sense designating by it digraphs with weighted 
adjacencies too. 

3 Circuit- and tree-like structures 

The matrix (graph) spectral analysis can be carried out using some form of character- 
istic polynomial of matrix itself or some special matrixes (graphs) constructed from it. 
Thus it is valid known "theorem on coefficients for digraphs" |@]. 
Theorem. Let 

det(AI - A) = + a^\^~^ + ■ ■ ■ + oat 
he characteristic polynomial of arbitrary digraph A with weighted adjacencies aij . Then 

a,= ^(-l)^Wvrz. , t = l,2,---,N , (1) 

Led 
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where Li is the set of all linear directed subgraphs L of graph A with exactly i vertices; 
p{L) means the number of components (dicircuits) of L. 

Coefficients in (|I]) are expressed in "circuit" structure of A, and this theorem is not 
more than rephrasing from the standard determinant notation |AI — A| in the form of 
matrix elements products sum with sign determined by substitution even (the number 
of dicircuits) into graph terms. 

In terms of "tree" structure it is known the characteristic polynomial expression not 
for matrix G itself but for its Kirchhoff matrix C (or conductivity matrix) determined 
hke 

C=D-G , 

where D is weighted powers matrix 

N N N 

D = diag(E gij. E 92j. " " " , E 9Nj) ■ 
Corresponding expression has a form : 



N 



det(AI-C) = (-l)^EM)'=[ E 

= n h = n ■ 

{i,j)£AF ii,j)eAF 

Here the set of forests containing directly k trees and being factors of G is designated 
by J^'^{G). Note that as the sum of elements along every line of C is equal to zero, so 
its determinant is equal to zero too and the sum in (refein) one can lead from k = 1. 
In the following in clear cases we omit indication on graph. 



4 The characteristic polynomial in tree-like struc- 
ture terms 

Knowing the characteristic polynomial expression of the admittance matrix C in terms 
of tree structure it is not hard to get analogous expression for characteristic polynomial 
of the matrix G itself. For this aim let us construct from G some new graph by 
the next rule. Let = U U {f} be the set of vertices A/" = {1, 2, ■ ■ ■ , A^} of G 
which is supplemented by some new vertex designated f. Let also add to G the vertex 

N 

t and lay on out of every vertex i G A/" the arc (i, f) with weight Qi] = — 9iji 

and remove all loops {i,i). We denote the obtained graph by G^ ( : VG"!" = N\ 
AG^ = {AG \ VJjyi^i)} U^(i, t) , weights of arcs are equal gij , i e Af , j E N"^). 

Corresponding generalized adjacencies matrix G^ of G^ has the following form 
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9lN 


911 ^ 


92N 


92] 


93N 


93] 





9n] 





J 



Gt 



9ni 9n2 9n3 

V 

Note that, if G = P — I, where P is probabihty matrix setting finite Markov's chain 
with kilhng, so the quantity is the probabihty of kiUing of the process if it is in state 
i. This has a sense of probabihty of outcoming to the bounder, so the additional vertex 
t could be interpreted in some sense as a bounder of the finite set A/" = {1, 2, ■ ■ ■ , A^}. 

It is easy to see, that (A^ + 1) x (A^ + 1) admittance matrix of graph has the 
form 



ct 





91] 


G 


92] 




9n] 


■■■ 






so, using (□) we get the chain of equations 



det(AI-G) = -det(AI + C^) 
A 



1 A^+l 



k=l 



N 

= EA'=[ E 

Note, that since the vertex f in graph G"^ has zero outdegree and hence it is a root in 
every forest F G JF'(G^), so the sets JF^+^(G'"'') and jFj^(G'^) coincide. Thus it is valid 

Theorem 1. Characteristic polynomial of an arbitrary N x N matrix G can be 
expressed in the form 



N 

det(AI-G) = EA'[E 



k=0 FeT'^ 



^F= n 



(3) 



where J^^ 



Under A = we obtain obvious 
Consequence. The determinant of N x N matrix G can be expressed in the form 



detG = (-1)^ E ^F 



(4) 
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Let 71 he a. subset of J\f. Designate by 07^.7^ matrix obtained from G by striking 
out the columns and lines with numbers i E TZ. So, 67^7^ is a diagonal minor of G of 
{N — |7?.|)-th order. The corresponding to it digraph we denote by G-jm. 

Consequence of consequence. The determinant of minor G-jm can be expressed 
in the form 

detG7^7^ = (-1)^-1^1 E , (5) 

Proof. By formula @ 

det Gnn = (-1)^-1^1 ^ . 

Let us keep in Gt^t^ the same numeration of elements as it is in matrix G. So the 
elements {gnn)ii of the corresponding matrix GIih (weights of arcs (i, I) of graph G\^-j^ 
are equal: 

{9nn)ii = 9ii , ^, / G A/" \ 7^ , 

{gimh = /} n 7^ ^ 0, 

{gnn)ii = E gim ,i & J^\7i . 
m&{-\}un 

So productivity np of any forest F G ^1(^7^7^) represents the productivity TVp' of 
some forest of the set jF{|}u7j(G"l'). The sum of productivities np along all forests 
F e ^t(G7j7^) exhausts the set J-'{^yuTi{G'^), which proves (^. 



5 Formulas for components of eigenvectors 

If it is known the eigenvalue A of G, so to calculate the components Vm of corresponding 
eigenvector v it is necessary to decide the standard system (AI — G)v = . Let for 
certainty the n-th component of eigenvector v be not equal to zero. Without loss of 
generality one can accept it be equal to one: Vn = 1. Then we have for the rest of 
components following the Kramer rule 

_ det A;,„,(A) 
detA„.(A) ' 

where 



A„„(A) — AI — Gnn 5 

G„n is algebraic adjunct of the element gnn of matrix G, and AJ^^(A) is a matrix 
obtained from A„„(A) by substitution of the m-th column of AI — G by n-th one with 
negative sign. 

The expression for A„„(A) it is easy to get using already obtained formula con- 
cerning diagonal minors Gt^t^ of matrix G. Since the coefficient at A'^' at the expression 
of characteristic polynomial of algebraic adjunct Gnn is itself a sum of determinants of 
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diagonal minors of matrix G of (A^ — /c — l)-th order with sign (—1)^ ^ and not 
including n-th column and n-th line so one can write 



N~l 



det(AI - G„,„) = E ^'i E det(-G^^) 

fc=o neAf, nen 
\n\=k+i 



N-l 

N-l 
k=Q 



-1 



.N-k-1 



E 



E det Gnn 

\n\=k+i 



E 



Hp 



\'R\=k+l 



The sum at square brackets in the last expression obviously is equal to the sum of 
productivities np of all forests belonging to the unification 

U ^{t}u7^(Gt) , 
7^eA/', n&n 

\-R\=k+l 

which in its turn is equal to J-'^^^^y{G''). So, the denominator at (||) has the form 



det A„„(A) = det(AI - G, 



^k 



{t,«} 



N-l 
k=0 



E 



■ Fp 



(7) 



It is more difficult to obtain the formula for numerator at (P) having the form 



Af-l 



det a; JA) = 5: 



fc=0 



E 



where by J-'^^^y{m ■ n) the subset of J^^:^^^{G^) having m ■ ri-walk is denoted. The 
summing over k one can lead till N — 2, because the set ^|^~|(m ■ n) is empty. 

It is sufficiently to establish (||) under A = 0, because the coefficients at powers of A 
are the diagonal minors of the same matrix and are utterly of analogical to each other 
form. 

Let A'^fn = ^nm(O)- necessary to prove that 



det = J2 t^f 



(9) 



In passing let us remark that using (j^) it is easy to obtain the expression for algebraic 
adjunct Gnm of the element gnm- Actually, matrices — A^^ and Gnm differ only in such 
a way. In — A^^ (if for example n < m) n-th column of matrix Gnm is at the (m — l)-th 
place and is of opposite sign. Hence 
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To prove let us exchange in G places of n-th and m-th columns. Obtained 
auxiliary matrix we denote by H. Nondiagonal elements hij of H and corresponding 
expanded matrix are equal 

i) hij = Qij , i ^n,m , 

ii) dim 1 him Qin 

It is easy to see that algebraic adjunct H„„ of the element of H differs from 
matrix (— A^^) by only the sign of m-th column. So taking into account 

detA;^ = -(-l)^-MetH„„ = - ^ . (10) 

Now in (|10D it is necessary to cross from the sum of productivities ttj? of factor forests 
of W to the sum of productivities vr^ of corresponding subgraphs of graph Notice 
that in W weight of arc (m, n) is equal to minus sum of weights of arcs (m, i) going 
out of the vertex m in graph G^: 

hmn 9mm ^ ^ Qmi ■ (-^-^) 

That is why we divide the set of forests in to two nonintersecting subsets: the 

set ^{i^njl-f^^; n)) of forests including the arc (m, n) and the set J^{^^n}{H^', 1 {m, n)) 
of forests without this arc. Then 

E ^F= E E ^F. (12) 

The set of arcs AF of the forest F G jF|| „}(i7^; (m, ra)) is representable in a form 
= {m,n) U , where P is some graph belonging to the set J-'{^^n,m}{H'^)- By it 
as hij = gij, j m,n the sets jFj^ „ „^}(if^) and J-'^-^^n,m}{G^) coincide. Hence, taking 
into account (|TT]) 



E T^F = - J2 9mi E 

Fe.F{t,n}{^^t;(m,n)) ieArt\{m} LFG.F{t,„,™} (Gt) . 

Note that any forest F G ^{t,n,m}(G'^) consists of exactly of three trees T^, Tn and Tm 
with roots correspondingly f, n and m. Therefore the addition to such forest the arc 
(m, i) with weight gmi leads either to graph belonging to J-'^^^n}{G^) (if ^ £ VT| U VT„), 
or (if i G VTm) to graph, in which outdegree of every vertex is equal to one with the 
exception of f and n: d'^{j) = 1 j ],n , (i+(t) = d^{n) = 0, and there is exactly 
one circuit in graph and this circuit contains the vertex m. We denote the set of such 
graphs by Om{G^)- So 

E 7iF = - E Tip- J2 ^F- (13) 
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Let us now consider the last sum at the right part of (|T2D. The productivity np of 
any forest F G „ ] (m, n)) by force of i) ii) one can represent as up = np, 
where E is some subgraph of G\ to obtain which from the forest F one must exchange 
the arcs ending in m into arcs ending in n and back. The forest F G ] (m, n)) 

consists of two trees T| and T„ with the roots f and n correspondingly. By this if 
m G VT-j- so the graph E is still a forest and besides the sequence of arcs from m does 
not lead to n (but leads to f). In the case of m G T„ in F such an exchanging of 
arcs results in E be graph containing the only circuit and the vertex m belongs to this 
circuit, i.e. E G (9m (C^). Thus 

E E 7rp+ . (14) 

Here we denote by jFj-j- „}((?''"; ]m ■ ra) the subset of J-'{-^^n}{G'^) not containing m ■ n-way. 
Combining together (^),(|1^),(|13]) and (|1^ we get, that products np along graphs 
containing circuits reduce and 

det = - E ^p= E ' 

which is the same as (^. Thus it is proved 

Theorem 2. Let A be the eigenvalue of order 1 of matrix G, v — corresponding 
eigenvector and its n-th component is not equal to zero, then accurate to constant 
factor the components of v are following 



fc=0 


E T^F 


Af-l 










E T^F 




A:=0 




^ppfk 
^■^{n.t} 





where fjl'^ ' ''^) ^■^ « subset of set consisting of forests having m ■ n-way. 

Remark. To obtain the expression for m-th component of eigenvector of the trans- 
form matrix it is necessary only to exchange indices m, n with each other at the 
right part of (|T5]). 



6 Example 

Let us consider an easy example demonstrating the application of singless form tech- 
nique. Let M = P — I, where P is a probability matrix, setting the finite Markov's 
chain with killing (that is the sum of elements along any line may be less than one, 
corresponding residual is just a kilhng probability), and also let the transition prob- 

N N 

abilities Mij = Py, i ^ j, and the killing ones = — J2 ^ij = 1 ~ E Pij be 

exponentially small. Namely, let = 3 and Mij = rriije'^^^^'^ and ^12=^13 = 4, 
Vi^ = 5, V21 = 3, V23 = 2, V2t = 5, V32 = 1, Vst = 4, V^si = 3; e is a small parameter. 



8 



Thus if instead of diagonal elements Ma we use their expression through the quantities 

N 

Mi^ {Ma = -Mil - J2 Mij) the matrix M gets the form 



-Mi2 - Mi3 - Mit 
M = I M21 
M31 



M12 

-M21 - M23 - M2t 
M32 



Mi3 

M23 

-M32 - M3t - M31 



By virtue of lack of sign and non-negativity of transition (nondiagonal elements of 
M) and killing (M^) probabilities that are used at tree-like structure formulas, it is 
ought to keep at the asymptotic of characteristic polynomial coefficients (|^) only terms 
reaching the maximal order on small parameter. Thus the eigenvalues asymptotic one 
can extract from the equation 

+ A^m32e~^/^ + Amgsmaie"^/^ + m32m2imi|e'^°/^ = . 
The exponential orders Vk = — lim e In of the coefficients ak of characteristic polyno- 
mial J2 A'^afc satisfy convex nonequalities system analogical to p: Vk-i — Vk > Vk — Vk+i, 
and the eigenvalues asymptotic one can look in the form A|. ~ e ? 0]- Substi- 
tuting model A in such a form we get the following eigenvalues asymptotics: 



Ai ~ — mi+e 



-5/e 



A2 ~ - m2ie 



, A. 



£^0 



"^326 



-1/e 



The asymptotic of eigenvectors of the matrix M and the transform matrix M* one 
can find from (^). Denoting by C (by C) the matrix, composed of ultimate values of 
vector- lines of M (vector-columns of M*) one gets 





(1 


1 






(1 -1 


c = 





1 


:) 




1 




[0 









i^O 







CC = I 



Notice also that the matrix 



M 



— mije 
m2ie^ 






—m2ie^~ 

1 



m32e c -^326 

demonstrates the same in leading order eigenvalues asymptotic and the same ultimate 
values of eigenvectors' components. 

Although matrices M and M outwardly are rather different, however the solutions 
of the evolution equation 



dp 



Lp 



with guiding matrix L equal to M* or to M* are the same under exponentially large 
times t = re^/^, V > 0, and in "slow" time" r with exponential scale V the corre- 
sponding evolution Pvi'^) = limp(re^/^) has the form 
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(r) = C' limdiag |e 



,AlTOxp(^) A2rexp(-5i^) Airexpi 



^"^^1 Cpv(O) . 



5 ) 



Under exponential scale V = 1 part of the initial distribution concentrated at the 
third state crosses to the second one. Such crossing rapidness at "slow" time r is 
determined by the quantity A3 = —77132- Under exponential scale = 3 if r — cx) entire 
distribution turns out to be at the first state, corresponding rapidness is determined 
by the quantity A2 = —m2i. Finally at the scale V = 5 killing of the process takes 
place and it is governed by the quantity Ai = —m^. Note, that in spite of the 
killing probability M^j = m^^e'^^^ at state 3 is greater than the killing probability 
Mij = mi-j-e"^/'^ at state 1, but the part of distribution concentrated at state 3 crosses 
to the second state (and then to state 1) before kiUing at state 3 could take place. 

Considered example shows from one side small matrix elements are not negligible 
comparatively with the large ones. From another side nevertheless one can neglect 
by some of the elements, may be not small ( here these are elements M12, M2-\, M23, 
M31, Msj). Which of elements do not affect on the spectrum and on the evolution is 
determined by signless formulas (||) and ([T5|) . 



vectors (p!5[) at the situation of arbitrary matrices give few substantial for practical 
computation in comparison with usual methods. Moreover, the formula (|l]) being only 
reformulation of standard characteristic polynomial expression is more preferable than 
@, because the volume of calculations needed for determination of coefficients at pow- 
ers of A using (JID is less than using (H). The reason is that G has less subgraphs being 
circuits than G"^ has forests. In such case signless formulas obtained are not of interest 
for practical calculations. 

However, the situation strongly changes in the case of large matrices with non- 
negative elements which coefficients depend on small (large) parameter and may possess 
on it different orders (such matrices, as it was noted before, often appear in singular 
and stochastic problems). Usual methods in such a situation are practically unsuitable 
due to large amount of calculations and low accuracy. On the contrary, formulas (|^) 
and ([T5| ) due to sign absence at sums become extremely effective and allow to keep an 
eye on only higher order in asymptotic. This circumstance permits not only to carry 
out practical calculations, but also (which is highly essentially from theoretical point oh 
view) to determine singular limits and analyze structure spectrum properties. Here it 
is necessary to note that though for matrix G with non-negative elements introducing 

N 

quantities = — J2 9ij are not positive however at the expense of spectrum shift 

it is always possible to make them non-negative ones. The problem of picking out 
at characteristic polynomial terms having highly order is reducing to determination 
forests F with maximal productivity irp and is a separate difficult one. This problem 
is solved, the effective technique for determination of extreme forests is elaborated . 
Analogical technique can be used for analysis of large very sparse matrices. 



7 Discussion 



Obtained expressions 




components of eigen- 
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